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ORIGINAL ARTICLE 

Scale Reduction of a Systems Coagulation 
Model With an Application to Modeling 
Pharmacokinetic-Pharmacodynamic Data 

A Gulati 1 , GK Isbister 23 and SB Duffull 1 

Bridging systems biology and pharmacokinetics-pharmacodynamics has resulted in models that are highly complex and 
complicated. They usually contain large numbers of states and parameters and describe multiple input-output relationships. 
Based on any given data set relating to a specific input-output process, it is possible that some states of the system are either 
less important or have no influence at all. In this study, we explore a simplification of a systems pharmacology model of the 
coagulation network for use in describing the time course of fibrinogen recovery after a brown snake bite. The technique of 
proper lumping is used to simplify the 62-state systems model to a 5-state model that describes the brown snake venom- 
fibrinogen relationship while maintaining an appropriate mechanistic relationship. The simplified 5-state model explains the 
observed decline and recovery in fibrinogen concentrations well. The techniques used in this study can be applied to other 
multiscale models. 
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Venom-induced consumption coagulopathy (VICC) is one of 
the most important snake envenoming syndromes. It is char- 
acterized by activation of the coagulation pathway by proco- 
agulant toxins in the venom, resulting in the consumption of 
clotting factors such as fibrinogen, factor V, VIII, and X. 1-4 A 
comprehensive systems pharmacology model of the coagu- 
lation network has been shown to describe the time course 
of changes in coagulation factors in response to various 
anticoagulants 5 and Australian elapid envenoming. 6 Brown 
snakes are the most common cause of serious snakebite in 
Australia, 7 and their venom contains a potent prothrombin 
activator. Brown snake envenomation results in VICC with 
depletion of fibrinogen, factor V and factor VIII. 1 ' 6 The time 
course of changes in fibrinogen, from any cause, is poorly 
understood, 8-12 making this an ideal setting to investigate the 
hypofibrinogenaemic states. 

An established mathematical model of the coagulation 
network 6 is described by 62 ordinary differential equa- 
tions (ODEs) and 178 parameters. See Supplementary 
Table S1 for the numbered states of the model with their 
corresponding species and initial conditions. ODEs for the 
model are also provided in the Supplementary Data (origi- 
nal_model.m). The values of the parameters in the model 
can be obtained from our previous publications. 5 ' 613 This 
model consists of multiple inputs, outputs, and modifiers of 
these input-output relationships. Based on any given set of 
available data on an input-output relationship, it is possible 
that while certain states and parameters in the model are 
relevant to the relationship, others are either less important 
or may have no influence at all. It is important to differen- 
tiate between parameters and mechanisms that would be 
informed by a new data set and those that would not. Elimi- 
nating the parameters and reactions that are not informed by 



the data would result in a simpler version of the model that is 
specific to that particular input-output relationship. In doing 
this, it is important to retain the mechanistic structure even 
for simpler models. 

Model order reduction methods for reducing the num- 
ber of states and parameters of dynamical systems that 
are defined by ODEs currently exist. 14 Lumping is a model 
order reduction method where the original states of the 
model are lumped or merged to a reduced number of 
pseudo-states resulting in a reduced number of equations 
and parameters but with effectively the same or similar 
input-output behavior. Proper lumping is a special case 
of lumping where each of the original states contributes 
to only one of the pseudo-states of the reduced system, 
thereby forming groups that retain a clear physical interpre- 
tation. Note there may be more than one state lumped into 
a pseudo-state, but no state will be lumped into more than 
one pseudo-states. The technique of proper lumping has 
been used to lump physiologically based pharmacokinetic 
models. 15-18 An application of proper lumping based on 
searching through the possible combinations of lumpable 
states has been described for NF-kB signaling pathways to 
demonstrate its use and performance. 19 

Even though the number of parameters in the blood coagu- 
lation model can be reduced using a model reduction method 
like proper lumping, it is important that the model is iden- 
tifiable to get precise estimates of the parameters. Various 
methods for the assessment of identifiability of linear and 
nonlinear models have been proposed. 2021 An informal, uni- 
fied approach using an information theoretic framework has 
recently been developed for simultaneous assessment of 
structural and deterministic identifiability for fixed- and mixed- 
effects pharmacokinetic-pharmacodynamic (PKPD) models. 
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The approach has been applied to Bateman and Dost 22 mod- 
els and to a combined parent-metabolite model describing 
the pharmacokinetics of ivabradine and its N-desmethylated 
metabolite, S-18982. 23 

Aims 

This work aimed to explore a simplification of a blood coag- 
ulation systems pharmacology model for use in modeling 
PKPD data. Specific aims included the following: 

1. To create a simplified model for exploring fibrinogen 
recovery after envenomation that mechanistically aligns 
with the coagulation systems pharmacology model 

2. To extract the simplified model for estimation purposes 

3. To assess structural identifiability of the simplified 
model based on the inputs and outputs available in the 
data set 

4. To develop a population PKPD model for fibrinogen con- 
centration-time data based on the mechanisms appar- 
ent in the simplified model 



RESULTS 

Simplification of the blood coagulation systems pharma- 
cology model 

The original 62-state blood coagulation systems pharmacol- 
ogy model was simplified to a 5-state model (Figure 1) that 
described the brown snake venom-fibrinogen input-output 
relationship. An in silico brown snakebite followed by an in 
silico antivenom administration at 4h resulted in a similar 
consumption-recovery profile for fibrinogen using the 5-state 
and 62-state models. During this process, some parameter(s) 
were adjusted heuristically by modifying their values manu- 
ally to refine the predictions of the different structure of 
the simplified model. The values of the parameter(s) were 
adjusted (manually) based on how well the predictions from 
the simplified model compared with the original model. 

Figure 2a compares the model predictions from the 
original 62-state model with the observed fibrinogen con- 
centrations post-brown snakebite; Figure 2b compares the 
model predictions from the lumped 5-state model with the 
observed fibrinogen concentrations post-brown snakebite; 
Figure 2c compares the model predictions from the 62-state 
and 5-state models; and Figure 2d shows the percent rela- 
tive differences in the predictions from the original 62-state 
and lumped 5-state models. It was observed that out of the 
five lumped states formed, fibrinogen and the brown snake 
venom absorption and plasma states were the ones that 
were initially left unlumped; the fourth and fifth states being 
factor I la and a single lumped state that was a group of the 
remaining 58 states. Table 1 summarizes the five lumped 
states with their initial conditions and their corresponding 
states of the original model that were lumped. Lumping the 
states further worsened the fibrinogen profile significantly. 

Extraction of the simplified model 

ODEs for the unlumped states in the 5-state lumped model 
were explicitly written by eliminating the reactions from 
the original ODEs that were lumped based on the proper 



lumping technique. Because state 5 was a lumped state 
that was formed as a result of merging of various states 
from the original model, its ODE had to be explicitly written 
as if it had been an unlumped state. The clotting factor that 
was most relevant to the brown snake venom-fibrinogen 
relationship represented the fifth state. With this in mind, 
the fifth state was therefore assumed to be dominated by 
factor II. The extraction of the ODEs of the simplified model 
from the ODEs of the original model resulted in reduction 
of the total number of parameters to 11 compared to 178 
in the original model. Because the ODEs were extracted 
from the ODEs of the original model, this version of the 
model was referred as the "extracted model." The following 
are the ODEs of the extracted model with their respective 
initial conditions: Eq. 1 , ODE for state 1 — venom absorption 
state; Eq. 2, ODE for state 2 — venom plasma state; Eq. 3, 
ODE for state 3 — fibrinogen; Eq. 4, ODE for state 4 — factor 
I la; Eq. 5, ODE for state 5 — lumped state. 
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Note that in these equations, the form of the model reflects 
a time-varying system rather than a nonlinear Michaelis- 
Menten process. Hence, the values of v and K do not 

1 max m 

necessarily reflect their implicit meaning in standard enzyme 
kinetics (see Wajima etal. 5 ). 

Identifiability of the simplified model 

Assessment of identifiability of the extracted model using 
Population OPTimal (POPT) design software found that fixing 
the two parameters, pLumpeotate anc | ^umpeotat^ resulted in both 
the predefined conditions for identifiability being met. Thus, 9 
parameters out of the total 1 1 parameters in the extracted 
model were identifiable. The remaining two unidentifiable 
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parameters, pLum P8 d_state anc | ^umpecstat^ arose f rom the lumped 
state and were required to be fixed at 147 nmol/l/h and 0.01/h 
, respectively. These values were assumed to be equal to the 
values of p" and a" in the original model. 



Venom dose 



Venom absorption state 


r, venom 
"abs 




Venom plasma state 


^venom 





-Lumped_state 




Lumped state 



^Lumped 



Figure 1 Schematic representation of the lumped 5-state 
coagulation model describing the brown snake venom-fibrinogen 
relationship. Dotted lines represent production or catabolic 
processes, dashed lines represent activation processes, and solid 
lines represent mass balance of the states of interest. 



Modeling the fibrinogen concentration-time data using 
the simplified model 

The decline and eventual recovery of fibrinogen after brown 
snake envenomation was described by the 5-state model. 
Testing additive, exponential, and combined error models 
for residual unexplained variability showed that the com- 
bined error model described the data the best. Compari- 
son of the models with between-subject variabilities on 
one structural parameter or more showed that thejriodel 
with between-subject variabilities on d" a , v ^ ped - sta,e >lla , 
^Lum P ed_siate >na ^ ana , ^ provided the best 

"' ' ' venom ~ 

description of the data. The parameter estimates for the final 
model are shown in Table 2. Individual fits obtained from 
the final model are shown in Supplementary Figures S1- 
S3. It can be seen that the model was able to capture 
the trend in the data for all individuals. Figure 3 shows 
the corresponding visual predictive check using the param- 
eter estimates from the final model. The visual predictive 
check showed that the model explained the observed 
data well. In comparison to the 5-state model, the original 
62-state model was only used to predict into the full data 
set, and this was not a model-fitting process but rather an 
overlay of the model predictions over the data. There are, 
therefore, no individual fits from the original model that can 
be shown. 




Figure 2 Comparison of original 62-state and lumped 5-state models, (a) Predictions from the original 62-state model (red line) overlaid 
on the patient data from the Australian Snakebite Project showing decline and recovery of fibrinogen after brown snake envenomation. (b) 
Predictions from the lumped 5-state model (green line) overlaid on the patient data from the Australian Snakebite Project showing decline 
and recovery of fibrinogen after brown snake envenomation. (c) Comparison of the predictions from the original 62-state model (red line) 
and lumped 5-state model (green line) showing decline and recovery of fibrinogen after brown snake envenomation. (d) Percent relative 
differences between predictions from the original 62-state and lumped 5-state coagulation models for brown snake envenomation affecting 
fibrinogen concentrations. 
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Table 1 Numbered states of the lumped 5-state coagulation model with their initial conditions and the corresponding states of the original 62-state model 



State number in 
the lumped model 



States 



Initial conditions 
(nmol/l) 



Corresponding states of the original model 



1 L Venom absorption state C 1L (0) = 0.0075 

2L Venom plasma state C 2L (0) = 0 

3L Fibrinogen C 3L (0) = 8,900 

4L lla C 4L (0) = 0 

5L Lumped state C 5L (0) = 7,900 



28 
62 
14 
7 

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 29, 30, 
31 , 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 , 42, 43, 44, 45, 46, 47, 48, 49, 50, 51 , 52, 53, 54, 55, 
56, 57, 58, 59, 60, 61 



See Supplementary Table S1 for the numbered states of the original 62-state model with their corresponding species and initial conditions. 



DISCUSSION 

This study explores the simplification of a complex and com- 
plicated systems pharmacology model and its application to 
modeling PKPD data. A comprehensive model of the blood 
coagulation network was simplified and used to model data 
from a large study of snake envenomations. The simplified 
version of the model was shown to describe the time course 
of changes in fibrinogen concentrations following brown 
snake envenoming, which was comparable to that with the 
original version of the model. 

Use of proper lumping as a model order reduction tech- 
nique in this study showed that a 62-state blood coagulation 
network model could be significantly reduced to a 5-state 
model that described the fibrinogen consumption and 
recovery after envenomation from a brown snake compa- 
rable to the original model. Brown snake venom-fibrinogen 
serves as an example of an input-output relationship to 
show the application of proper lumping to reduce the num- 
ber of states in a systems pharmacology model. The tech- 
nique, however, can be applied to any other input-output 
relationship. 

The ODEs of the reduced system had to be written explic- 
itly, and the resulting extracted model only describes the 
brown snake venom-fibrinogen relationship. Lumped mod- 
els are specific to the input-output being studied and other 
lumped models will be required from the same systems model 
for different input-output relationships. Use of an Information 
Theoretic Approach to assess the simplified model resulted 
in the list of parameters that were identifiable and thus could 
be estimated precisely. The structural identifiability method 
used in this work was based on that proposed by Shivva 
et al. 23 This method is a local identifiability method based on 
an information theoretical framework. This could be termed 
as a "pragmatic" identifiability solution. 

In this study, the half-life of fibrinogen was estimated to be 
1 .5 days. This compares to a value of 1 day from the origi- 
nal 62-state model (as per Wajima et al.). 5 Also, this value 
is similar to another study of taipan bites in which the half- 
life was found to be 1 day. 24 Others have found longer val- 
ues of half-life between 4 and 5 days, 25 and it may need to 
be explored whether there are other components of snake 
venom that may eliminate fibrinogen via other mechanisms 
or there are other processes that get stimulated following 
snake envenoming. The half-life of brown snake venom was 
estimated as equal to 55min. There are no previous studies 
of brown snake venom in humans that have investigated the 



Table 2 Parameter estimates for the final model for fibrinogen 



Parameter 




Estimate 
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0.0174 


9.4 




pfs (nmol/l-h- 1 ) 




216 


6.8 




"mil * F ( nmo| / | - n_ 


') 


8.88 


4.5 




lla 

^Fg >F ( nm0 |/|) 




22.9 


5.7 








0.744 


0.87 








0.854 


24.2 




BSV _ d"" (%) 




0.00500 


10.0 


75.0 


DOW Lumped state — 
DO V _ V 


>"> (%) 


2.10 


23.0 


43.0 


BSV ^-Lumped_state — 


'""(%) 


0.990 


24.5 


61.0 


BSV_d Fg ( % ) 




0.530 


17.8 


46.0 


BSV _ F (%) 




1.40 
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2 venom 

a 
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0.390 


3.0 
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2 

add 




0.00100 


1.3 


9.56 



half-life of the venom. However, in other snake species like 
red-bellied black snake, the half-life has been reported to be 
-8-1 2 h (G.K. Isbister, personal communication). This was 
only able to be determined because the antivenom was not 
given in many cases of red-bellied black snakebites. In addi- 
tion, the half-life refers to the presence of venom as a whole. 
The half-life of brown snake venom estimated in this study 
refers in particular to the activity of the prothrombin activator 
in the venom that is Xa:Va like. However, it is possible that the 
half-life of the functional activity of the procoagulant toxin is 
different to the half-life measured as the presence of whole 
venom in serum. 

The use of lumping principles has been illustrated previ- 
ously for systems models in the area of PKPD. These works 
have concentrated on physiologically based pharmacokinetic 
models, and only limited attention has been placed specifi- 
cally on systems models describing pharmacological effects. 
The studies on physiologically based pharmacokinetic mod- 
els characterized the pharmacokinetics of a homologous 
series of barbiturates in rat 17 and for compounds like 1 ,3-buta- 
diene 15 that are largely used in the production of plastics and 
synthetic rubber. Proper lumping principles have also been 
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5 



4 




Time post-bite (h) 

Figure 3 Visual predictive check for the final model for the decline 
in fibrinogen concentrations and its recovery following a brown 
snake envenomation. Red lines represent the percentiles of the 
fibrinogen data and black lines represent the model predictions. 50th 
percentile is represented by the solid line compared to 1 0th and 90th 
percentiles, which are shown by the dashed lines. 

applied to a systems biology model describing the signaling 
pathways of NF-kB. 19 These studies, however, only discuss 
the illustration of the lumping process. The work carried out 
in our study appears to be the first time that a systems phar- 
macology model has been simplified to model PKPD data. 

It might be theoretically possible to perform an identifiabil- 
ity analysis on the original model and then fix unidentifiable 
parameters. However, this model (i) would be cumbersome to 
use and computationally intensive to model and (ii) would not 
be able to be described simply in any report or future work. 
The theoretical benefit of the methods presented here is the 
observation that a model can be "extracted" from the original 
model to accommodate potentially any input-output relation- 
ship that is of interest. In addition, the use of mechanistic ele- 
ments in models provides for greater predictive performance 
to settings outside of which the model was created. Using 
the simplification techniques discussed in this study, any 
multiscale mechanistic model that includes the necessary 
input-output relationship of interest can be used as a start- 
ing point for model extraction and then building. The result- 
ing simplified model preserves the mechanism as seen in 
the original systems model. This allows for rapid input-output 
model building that effectively eliminates lengthy structural 
model selection processes. However, it should be noted that 
model order reduction techniques are locally dependent on 
the parameter space. Care should be taken to ensure that the 
performance of the final reduced model meets the intended 
needs. 

To conclude, a simplified version of a multiscale model was 
created for exploring fibrinogen recovery after snake enven- 
omation. The simplified model mechanistically aligned with the 
coagulation systems pharmacology model and was extracted 
for estimation purposes. The simplification process is mul- 
tistaged and could benefit from some level of automation. 
A population PKPD model for fibrinogen concentration-time 



data based on the mechanisms apparent in the simplified 
model was developed and was able to describe the recovery 
of fibrinogen following brown snake envenomation without 
the need for further structural model development. 

METHODS 

Simplification of the blood coagulation systems pharma- 
cology model 

The technique of proper lumping based on a previously pub- 
lished method by Dokoumetzidis and Aarons 19 was used 
to simplify the blood coagulation systems pharmacology 
model. The work of Dokoumetzidis and Aarons consisted 
of two components (i) a formal method for lumping and (ii) 
an automated algorithm for lumping a systems model. The 
formal method of lumping demonstrated the use of a lump- 
ing matrix to simplify a model with X numbers of states to a 
model with V numbers of states (where Y< X). In our work, 
we used their formal method of using the lumping matrix, 
but instead of using their automated algorithm, we manu- 
ally undertook the process to control more easily what states 
were lumped. 

The formal method of using the lumping matrix to reduce 
the number of states of a model, as explained previously 
by Dokoumetzidis and Aarons, 19 is described briefly below. 
In our work, the resulting simplified model was termed the 
"lumped" model. The 62-state systems pharmacology model 
was referred to as the "original" model. Supplementary 
Table S1 summarizes the 62 states in the original model 
with their corresponding species and initial conditions. 
ODEs for the original model are also provided in the 
Supplementary Data (original_model.m).The values of the 
parameters in the model can be obtained from our previous 
publications. 5 ' 613 

The 62-state model was written in terms of the concentra- 
tions of the species in their respective states given by vector y 



where y is the vector of states of length 62 and t is time, 
an independent variable. 

A lumping matrix M consisting of 0s and 1s and of dimen- 
sion n L xn (n L < n) was used to transform the vector of states 
y of dimension nx1 to a vector of lumped states y L of dimen- 
sion n L x1 : 

y, = My and therefore = M— (7) 
L dt off 

The number of rows in the lumping matrix M was equivalent 
to the number of final lumped states intended to be obtained 
as a result of proper lumping. The number of columns in M 
was equal to the number of states in the original model which 
was 62 in this study. As an example, if the original 62-state 
model was being lumped to a 35-state model, dimensions of 
y L would be 35 x 1 and that of y would be 62 x 1 . This would 
give dimensions of M equal to 35x62, which would be used 
for the transformation of y to y L . 

Inverse transformation of Eq. 6 is given by: 

V = M + y L (8) 
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Here, WT is the Moore-Penrose pseudoinverse of the 
lumping matrix M.Thus, given the 62-state coagulation model 
— = f (y) and the lumping matrix M, which transforms the 
vector of states y to a vector of lumped states y L such that y L 
= My, the aim was to determine f L (y L ), where _ ^ ^yj. 

The following steps were followed to determine f L (yj\* 



ofy 
dt 



f(y) 



M 



Y = M + y L 

dt y ' dt 



Mf(M + y L ) 



Thus, the main lumping formula that was used to simplify 
the 62-state coagulation model was (as per Li and Rabitz 26 ): 



f L (y L ) = Mf(M + y L ) 



0) 



For linear systems, the lumping matrix M can be used to 
transform the parameter vector of the original model (K) to 
the parameter vector of the lumped model (K L ). The lump- 
ing formula used gives the lumped system explicitly with the 
values of the parameters of the lumped system. The lumped 
model can then be used for estimation purposes. For time- 
varying linear or nonlinear ODE systems, the lumping formula 
does not provide adjusted parameter values for the lumped 
model. But it gives the initial conditions of the lumped states 
that are formed as a result of merging of various states. The 
lumping process uses the following relation to give the initial 
conditions of the lumped states: 



y 0 , =My 0 



(10) 



Here, y„ and y 0L are, respectively, the vectors of initial con- 
ditions for the original and lumped models. Thus, in the case 
of time-varying or nonlinear systems, the lumped model pro- 
duces outputs based on the reduced number of states and 
modified initial conditions but with the parameters from the 
original model. 

For the time-varying nonlinear coagulation model, the 
above process was used to simplify the coagulation model up 
until the point when the lumped model started to give diver- 
gent results to the original model using the parameter values 
of the original model. After that point, the lumping process 
was not able to simplify the coagulation model further using 
the same set of parameter values as the original model. Fur- 
ther simplification of the model was only possible heuristically 
by adjusting values of parameter(s) manually. The values of 
the parameter(s) were adjusted (manually) based on how 
well the predictions from the simplified model compared with 
the original model. 

Based on the input-output relationship considered in this 
study, states that were left unlumped were fibrinogen and 
the brown snake venom absorption and plasma compart- 
ments. For the remaining 59 lumpable states, each state 
was lumped with another arbitrarily chosen state, and the 
lumping matrix M was constructed accordingly. This search 
for the combination of states was conducted in a stepwise 
manner looking through a single series of lumping steps. 



If we intended to lump our original 62-state model to just 
52 states, there will be 10 26 possible combinations. This 
total number of combinations was based on the binomial 
coefficient: 



f62x(62-1)/2'| 
10 



(62x(62-1)/2)! 
10!((62x(62-1)/2)-10)l 



(11) 



Ultimately, an automated global search algorithm is desir- 
able that operates within a predefined set of constraints. But 
this is not yet available. As the number of possible combi- 
nations of states was extremely high in this study, all pos- 
sible combinations were not tested in this work. The lumping 
method used here can, therefore, be considered as a local 
solution. 

The initial conditions of the lumped states obtained from 
the lumping formula and the parameters of the original 
model were used to simulate the time course of fibrinogen 
post-brown snakebite. The simulated time courses were 
then compared among the lumped and the original mod- 
els. To simulate the time course of fibrinogen in a complete 
VICC scenario post-brown snakebite, a virtual brown snake 
venom dose of 0.0075 nmol/l, as in a previous study, 6 was 
used. All complete VICC patients were given antivenom in 
the observed snakebite data set, 1 and its administration in 
the model was set at 4h after the bite as the average found 
in the brown snake venom data. Simulations were carried out 
using MATLAB R201 1 a. The simulated fibrinogen concentra- 
tions from the lumped model were compared with those from 
the original model at each time point to assess for loss of 
predictive performance. The lumped model was considered 
as giving a similar output to the original unlumped version 
if the percentage relative difference was less than 20% for 
all time points. The value of 20% was chosen arbitrarily. For 
differences greater than 20%, the state was unlumped from 
the pseudo-state being considered at that point and other 
combinations assessed. This process was stopped when no 
further combination of states was possible that resulted in a 
similar time course of fibrinogen post-bite among original and 
lumped models. 

Extraction of the simplified model 

Because the aim of this study was to reduce the number 
of states as well as parameters, it was important to extract 
the ODEs of the lumped model to model the fibrinogen 
concentration-time data. ODEs of the lumped model were 
extracted from the ODEs of the original model by eliminat- 
ing the "unwanted" reactions that did not have any influence 
on the fibrinogen profile based on the proper lumping tech- 
nique. The ODEs of the lumped states (or pseudo-states) 
that were formed as a result of merging of various states 
from the original model had to be explicitly written as if they 
had been unlumped states. The input-output model with the 
reduced number of states and parameters was referred to as 
the "extracted model." 

Identif iability of the simplified model 

The structural identifiability of the extracted input-out- 
put model was assessed using an Information Theoretic 
Approach. 23 POPT design software was used to carry 
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out the structural identifiability analysis. A venom dose of 
0.0075 mg and a generic study design ( I ) for the sampling 
times (1, 2, 3, 4, 5, 6, 8, 10, 16, and 24 h) were assumed 
for the analysis. A predefined condition, as per Shivva et 
al., 23 had to be met for a model to be structurally identifi- 
able. According to the condition, for all values of design 
variables (£), where number of time points of observation 
is greater than the number of parameters and all design 
variables have a unique value, the log of the determinant 
of the Fisher information matrix should approach infinity 
as the associated noise, the log of the residual variance, 
approaches zero. If the predefined condition was not met, 
different sets of parameters, individually and in combina- 
tion with each other, were fixed and the model reassessed. 

Modeling the fibrinogen concentration-time data using 
the simplified model 

The data. The data from the Australian Snakebite Project con- 
sisted of timed clotting factor concentrations in a large cohort 
of patients with VICC resulting from Australian elapid enven- 
omation. There were 1 38 snakebite cases recruited to the Aus- 
tralian Snakebite Project from over 100 hospitals in Australia 
between January 2004 and May 2008. Each patient contributed 
only one case. Table 3 summarizes the patient demographics. 
The data used for this analysis consisted of timed fibrinogen 
concentrations in patients with complete VICC resulting from 
brown snake envenomation. All patients with complete VICC 



Table 3 Patient demographics 



Total number of patients bitten by brown snake 


73 


Number of patients with complete VICC 


61 


Percent of male patients 


74% 


Age of the patients in years (mean ± SD) 


41 ±18 


VICC, venom-induced consumption coagulopathy. 




Time post-bite (h) 

Figure 4 Patient data from the Australian Snakebite Project showing 
decline in fibrinogen concentrations after brown snake envenomation 
followed by recovery after antivenom administration. Each dashed 
line represents data from a single patient. Circles represent single 
data points. 



had been administered antivenom. The data for many patients 
were sparse, and uncertainty in the venom "dose" was a key 
feature of the data. The limit of quantification for fibrinogen was 
0.2g/I.The data are shown in Figure 4. 

Data analysis using a population approach. A full popu- 
lation approach was carried out to analyze the fibrinogen 
data using NONMEM version 7.2 (using first-order condi- 
tional estimation with interaction). 27 Beal's method M3 28 
was used for handling data that were below the limit of 
quantification. The extracted structurally identifiable simpli- 
fied coagulation model was used as the structural model 
to describe the fibrinogen concentration-time data. The 
compartmental structure of the model was not subject to 
further modeling. The dose of brown snake venom was 
set to 0.0075 mg for all patients. Additive, exponential, and 
combined error models were tested for the residual unex- 
plained variability. Residuals were assumed to be normally 
distributed with a mean of zero and variance a 2 . Between- 
subject variability on the structural parameters was mod- 
eled assuming an exponential structure and was included 
on the bioavailable fraction, where the fixed-effects value 
was fixed to 1 . 

Any parameters determined to be unidentifiable were fixed 
at the values in the original 62-state model and thus not esti- 
mated. Model comparison was based on the likelihood ratio 
test with a decrease in the objective function value of at least 
3.84 for one parameter difference between nested models to 
determine statistical significance. Additional model selection 
criteria included reduction of random residual unexplained 
variability and the stability of the model. A percentile visual 
predictive check to evaluate the final model was performed 
by simulating 1,000 replicates from the model and compar- 
ing the 10th, 50th, and 90th percentiles of the observed data 
and the prediction intervals derived from the simulated data 
graphically. Censoring of the data was used with the simu- 
lated data for the visual predictive check in the same way as 
the real data. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y Existing systems pharmacology models are 
used for simulation and in some cases estima- 
tion purposes. Because they are quite complex 
and cumbersome to work with, estimation is 
generally not considered as their primary ap- 
plication. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

y This study addressed how multiscale models 
can be simplified and then used as a mecha- 
nistic basis to create models to estimate param- 
eters that describe an input-output relationship 
of interest. 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y This study applies existing techniques. We de- 
scribe the use of a formal model order reduc- 
tion technique to create a mechanistic simplified 
model that can be applied directly to modeling 
data. This mechanistic systems model can be 
taken as a starting point, therefore avoiding the 
need for building a model empirically. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y The use of scale reduction techniques de- 
scribed in this study would encourage the use of 
mechanistic models in the field of clinical phar- 
macology, the further development of systems 
pharmacology models, and provides a tangible 
method of using these models in clinical drug 
development. 
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